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^. 1 Introduction 

en 

p^ The CIR process has been proposed by Cox, Ingersoll and Ross [T] to model the evolution of 

. . interest rates. It is defined as the unique solution of the following stochastic differential equation: 

"X dVt = ia-bVt)dt + ay^tdWt, Vo = v, (1.1) 

where a, a, v > and 6 G M (see [6] for the existence and uniqueness of the SDE). The density 
transition Vt is known and is given by (cf. [7J) 



«("'^' = 2^(^j ""-Ar^i^yA:^/^)' "■'' 



^ I would like to thank Prof Damien Lamberton for many useful discussions. 



where a = - — ^ — , i' = 2a/ a and lu is the modified Bessel function of order z^. We set 



ft 
It 



/ Vudu (1.3) 

Jo 



GuHsashvili and Stein |1] (cf Lemma 6.3 and Remark 8.1) show that It admits a distribution den- 
sity (see also Dufresne |3] for a similar result). A sharp asymptotic formulas for the distribution 
density of It is obtained in |3| (cf Theorem 2.4) as 

piiy) = At e-^^J'+^'v^ yf-t (l + O(y-i)) (1.4) 

where A, B and C are positive constants. The integral of V arises in both volatility and interest 
rate models. For example, the CIR process is used to model the instantaneous variance of 
underlying asset in the Heston model |5J , where the dynamics of the log-price of the underlying 
asset is given by the SDE 

dXt = -h't + VVt [pdW} + Vl - pHWf'^ (1.5) 

whith p g] — 1,1[. Heston's model is certainly the most popular among stochastic volatility 
models. Its popularity is mainly due to the availability of semi-explicit formulas for the prices 
of European options on the underlying asset using the fact that the Laplace transform of Xt is 
given explicitly. However, there is no explicit formulas for the distribution density of Xt. In the 
particular case where p = 0, a sharp asymptotic formulas for the distribution density is given in 
0] as 

pf'P=^{x) = Aie-^^^+^^V^x^^J (l + 0(x-3.)) (1.6) 

where ^i, A2 and A^ are positive constants. In general case, Dragulescu and Yakovenko |2j show 
that the density of Xt behaves as e~^+^ for x — )• -|-oo (and e~^~^ for x — )• —00). In particular, 
the moment generating function of Xt becomes infinite ('explodes') after the 'critical moment' 
p,* = p+. This moment explosion phenomenon has also been linked to the wing behavior of the 
implied volatility, via Lee's celebrated moment formula [8]. 

The purpose of the paper is to connect the tail behavior of some real-valued random variable 
to the behavior of the moment generating function near the 'critical moment'. Our main result 
is to show that if Z is real-value random variable for which /i 1 — > E e^ explodes at some critical 



moment //* and behaves in a certain way near p* (see Assumption 3.1), then one can obtain a 
sharp asymptotic formulas for the cumulative distribution of Z, where the "secondary-terms" 
are given in limsup statement. We show that this result can be applied to random- variables of 
type AiVt -|- \2lt as well as Xt. As an application of our results, sharp asymptotic formula for 



the implied volatility in the Heston model is obtained. 

This paper is organized as follows: In section 2 we give the moment generating function of 
an arbitrary combination of V and /. In section 3 we give a large deviation result linking the 
moment explosion and the tail behavior. The section 4 gives the tail behavior of X and the wing 
behavior of the implied volatility. We give some proof in the appendix. 

2 Moment generating function of a combination of V and / 

To compute the moment generating function of some combination of V and /, we use the addi- 



tivity property of V with respect to (a, v). Indeed, if we denote by V^'^'^the solution of ( 1.1 ) with 
initial condition v, then for any (fi, f2),(ai) ^2) G I^+j the process Y^i+^2,ai+a2 j-^g^g ^j^^ same 



low as (V^^^''^^ + V^^'""^) where Y'"''^ is the unique solution of (1.1) in which W is replaced by an 
independent Brownian motion. Indeed, by setting Vf^ = y^'^''^'^ ^ y^ = ^^2,02 ^^^ ^ _ y^i _j_ y2^ 
we have 

Vt = V1+V2+ I [(ai +a2)-hVs]ds + a I ( ^/v}dWs + ^^dWA 

= V1 + V2+ I [(ai + 02) - hVs] ds + a y^sdBs, 
Jo Jo 

where B is Brownian motion defined as 

dBs = , ^ ' dWs + , ^ dWs. 



It follows from that if we set, for Ai and A2 G M„ 

we have F^i+"2(t,ui + ^2) = F^ ^^{t,vi)Fll ^^{t,V2). So the function F^^^^^ takes the form 

The functions ip and -0 are characterized by the following result, whose proof can be found 
in the appendix. 



Theorem 2.1. For Ai, A2 G K, denote by V'Ai,A2 t^^ maximal solution of 




(2.1) 



defined over [0,t*^ ^ [ and denote by ip\^^\^{t) := L 'ip\-^^x.^{u)du. Then, for any T < t*^ ^ , w;e 
have 

Jg gAiV:^'''+A2/T _ ga¥'Ai,A2(T')+«^^Ai,A2{^) (2.2) 

The function F^ -^ is defined over [0, t^ _;^ [ s.t linii_i>i* F^ ^ {t^v) = +00. For Ai, A2 G 

M, denote by tl^^^^{fi) = t*Ai,/xA2- Denote also by fi*_{t) = -sup |^ < : tl^^^^{fj.) = tj and 

^+(i) = inf < /i > : t"^ X (/u) =t( (with sup$ = —00 and inf<l> = +00). So the moment 
generating function of AiVi+A2/t is defined over ]—fi*_{t),fi\{t)[ and if fi*_{t) < (resp /U^(t) < 0) 
then Hm^^_^;^(t) E e^^i^'+^^a-^* = +00 (resp hm^^^* (j) E e^^i^*+^^2Jt = _^^y 



Proposition 2.2. Lei i > and Ai,A2 G M so that max(Ai,A2) > 0. Then /i+(t) < +00. Fur- 
thermore, z 
near fJ.'^{t). 



thermore, 3w+ > s.i // 1 — > An (e e/^^i^*+/^^2 /q v„d«\ _ ( ,'^J^_^ - ^In ^, ^]-^_ )] is bounded 



3 Moment explosion and asymptotic behavior of the cumulative 
distribution 

In this section we present our main result based on the moment explosion of some random 
variable Z 

Assumption 3.1. There exist fi*,uj G M;^ and u £ R so that E e'^^ < +00, V /i G [0,^*[ and 
(in (E e^^) — A(/x)) is bounded near jjl* , where A(/x) = ,'^ + z^ln }_ ■ 

We next present our main result concerning the asymptotic behavior of the cumulative dis- 



tribution of a random variable Z for which Assumption 3.1 holds. 



Theorem 3.2. Let Assumption \3. 1\ hold. Then we have 






\n¥(Z> R)+ii*R-2yMR 
hmsup — -— G 



V 3 v 
2 ~ 4' 2 



(3.2) 



turtliermore, if iim/j_!>oo — ^ irrrsN = ^ . tlien 



In(iJ) 



V 3 

2 ~ 4 



(3.3) 



Corollary 3.3. For any t > and Ai,A2 G M so that max(Ai,A2) > 0, p^ip , (g.gp and (g.gp 
/lo/d for Z := A1I4 + ^2h- 



4 Asymptotic behavior of the stock price distribution in Heston 
model 

The analysis relies on the explicit calculation of the moment generating function of X. Let's 
define the quantity 

Fp{t)=EeP^', peR. (4.1) 

We have for all p G M 

F (t) = Ee-2fo Vsds+pp /o' VV^dW^ +P,/i^ /o* VVsdW^ 



EU 



I^Sl^p^ /„' Vsds+pp f^ VVsdW^ 



EeP- 



./T^f^VVsdw!- 



,2_r^(i_phjty^a^ 



(wr. 



s<t 



E 



ppj^, VV^dW^-^J^ Vsds J-^ Ji Vsds 



E^e'^i'o^'"^^ 



where we used the law of iterated conditional expectation and the fact that Vg is measurable 
with respect to W^. The last inequality is a consequence of Girsanov theorem, where under Q, 
the process V satisfies the stochastic differential equation 



dVt = {a-{b- p(jp)Vt) dt + a^JVtdWl 



(4.2) 



with Q— Brownian motion W^'^. We are the reduced to the calculation of the moment generating 
function of the time average of the CIR process V under Q. It follows that Fp{t) = e«¥'p(*)+''V'p(<)^ 
where (pp{t) = L 'ipp{s)ds and ipp is given by theorem 



2.1 



Theorem 4.1. For any t > 0, Assumption\3. 1\ holds for Xt as well as for —Xt. In particular, 



(3.1), (3.2) and (3.3) hold for Xt and {-Xt). 



Asymptotic behavior of the implied volatility We can also analyze the asymptotic be- 
haviour for the implied volatility. Recall that the implied volatility at = crt{k) of a call option 



St = e * with strike K = e , and maturity t is determined from the relation: 



E(e^'-eM = CBs{t,k,at) := N{di{t,k,at)) - e''N{d2{t,k,at)) 



(4.3) 



where 



di{t, k,a) 



^^ and d2{t,k,at) = di{t,k,at) — Via. 



a 



and N is the cdf of the normal law AA(0, 1). The asymptotic behavior of the implied volatility is 
linked to the tail behavior of the distribution of the stock price process. Here we get the following 
result 



Theorem 4.2. for any t > 0, we have 



lim 

k—^+oo 



at{±k)-(3fVk 



ki 



lim sup — 

k^+Qo /^g ln(/c) 



at{±k)-{pfVk + p^)) e 



0, V7> 



a 1 a 1 

^2 "4'a2 + 2 



(4.4) 
(4.5) 



And if (3.3) holds for ±Xt, then 

Vk 
lim 



fc-s>+oo f]f ln{k) 



at{±k)-{P^Vk + /3f) 



a 1 



(4.6) 



Whzthpt = ^ [^|^-^Jpl - l) , p.. 



'2 -^jr 



Ip^-1 Jp^ 



P: 



X ' ^3 - V2i 



/p±-l Jp± 



p+ = ^+ (t) and p_ = lJ._{t) + 1, where fM^ and iJ± are given by (E.4), (E.5), (E.l), (E.IO) 



and (E.ll). 



A Proof of Theorem 2.1 



Denote by 



F^ ^ (t,t,) = e°"^^l'^2(*)+^V'Ai,A2(t)_ 



The function F^ ^ is the maximal solution of 



Ai,A2 



2 fp,Tpa r)F^ 

-(*'«) = ?r^^^(^'^) + (« - hv)^^{t,v) + A2t;F? ,,(t,t;), < t < T, 



9t '^'"' 2 ■' 9^2 



dv 



[ F^x,MiO,v) = e'^^ 



Indeed, this is equivalent to ipx-^^x^{0) = 0, ^Ai,A2(0) = -^i ^-^d 

2 
Vai,A2(*) + ^V'li.Aa (*) = Y^'V'!i,A2(*) + (« - &l')V'Ai,A2(i) + >^2V 

= aV^Ai,A2(i) + ^'(yV'i,A2W-^V'(t) + A2). 
This is equivalent to ip'^_^ ^^{t) = ^Ai,A2(0 ^^^ 

i^'x.Mit) = y (v'!„A2(i) - 2^V'Ai,A2(t) + 2^ 



So for T < t*(Ai,A2), the process (e^^fo^^''''^'^F^^^^^{T-s,Vs''')) is a positive local 
martingale. It is therefore a super-martingale. In particular, for any stopping time r G 7q,t, we 
have 

It follows that 

sup ii; e'^2/J'K"'"rf«+«</'Ai,A2(T-T)+yT-"'"VAi,A2(T--r) ^ ga<^Ai,A2(T)+«V'Ai,A2(T) (A.l) 

We will show that there exists p > 1 such that 

sup ^eP(^2/,;K"'"d«+aipAi,A2{7^-^)+Vr'>Ai,A2(7^-^)) <+oO 

Note that ^Ai,A2 is increasing with respect to Ai and A2. This means that if Ai < Ai (resp 
A2 < A2), we have Vai,A2(*) < V'Ai,A2(*)' ^o^ '^^'^^y * < t*{^i^^2) (resp V^Ai,A2(i) < V'Ai,A2(0' 
Vt < t*(Ai, A2)). Indeed, if we set / = 'ij^x^^x^ " V'Ai,A2 and g = ^x^M ~ V'Ai,A2, we have 



/(O) = Ai-A2>0, 

fit) = y(^i,A2W-^l,A2-^/W 

= y/w(^Ai,A2W+V'Ai,A2-^ 



and 

9{0) 



9' it) = y (<,.,(*) -<A.-^5W + 2- 



2^ ..X <.A2-A2 



^2 



> y5(i)(^V'Ai,A2(0+V'Ai,A2 ^2 

It follows that / and g satisfy (/e^)'(t) = 0, with /(O) > and (5e^)'(t) > with ^(0) = 0, 
where A'{t) = ^ (Va^.a^C*) + V^Ai,A2 - f )■ Thus / > and 5 > 0. 
It follows that 

V'Ai.A.W - V'A.A.W = (Ai - Ai) e''*-i'o('^^i.-2+V'Ai^^2)(-)'^«, Vt < t*(Ai,A2) 

So for every T < t*^ ^ and for e > small enough, there exists AJ > Ai such that 

{l + €)^Px^Mit)<iJ^lM{t), Vi < T. 

We deduce that for any stopping time r € Tq^t, 

where a^ = (1 + e)a et A| = (1 + e)A2. Thus 



reTo^ 



T 



In conclusion, the local martingale (Ms := e^^o ^" dupa ^ ^p _ s^Vs'^)) is uniformly 
integrable; it is a true martingale. In particular, we have 

E eA2/o*\C-"d«+a^;,^,;,2(^-*)+^t'''>Ai,A2(T-t) ^ ^av^Ai .Aa (T)+V;"'>Ai.A2 (T) ^ y^ < 2._ 

Thus 



B Solution of (|2J^ 



To solve (2.1), we write 



6 , ,.. . „A2 /, , , bV 2A2fT2-62 



a^ 



G^ 



i^XiA'zit) 



+ 



so there are three situations : 
Case 2A2cr^ < 6^ : By setting 



ipx^Mi't) = V'Ai,A2(*) o and q = v^P^^2A^/ct^, 



the equation (2.1) becomes 



(7^ 



a 



V^l„A.W = v(^A„A.W-«^) 



which gives 



Then, 



It follows that. 



2a VV'Ai.AaCO -a '0Ai,A2(i)+a 



In 



V'Ai,A2(0 -Q 
V'Ai.AalO + a 



QO" i + In 



i^'x.Mit) 



Ai-4 + a 



O" 



V'Ai,A2(^) = -« ^g2^2..2,J^ , SI |Ai-^|>a, 



^ V'Ai,A2(0 = -« 'g2e2..^t_i ^ Sl|Ai-^|<a, 



where C 



Ai - ;^ - a 



Ai - T^ + a 



Thus 



V'A„A2(t) = ^-«g^, if |Al-;^|>a, 



V'Ai AaW = 4 - a "-^ 2r\ if |Ai - 41 < a. 



Case 2A2cr = b : The equation (2.1) can be written as 



a 



i^X.Mit) = vV-'Ai.AaW 



which gives 



Thus 



V'Ai.AaC*) 






2 



+ 



Case 2A20"^ > 6^ : We set /3 = ^2X20'^ — h'^ jo^. The equation (2.1) can be written as 

J2 



A.M^t) = - [^\^^^{t) ^ f) 



So 



^M,A.(^) 



(T 



Integrating both sides of this equation we obtain 



1 



oH 



arctan(V'Ai,A2 (*)//?) = "Y" + ^*^- 



Thus 



^Ai,A2(i) 



a^ 



+ /3tan ( /3 h arctan 



Aia^ - h 



/3a2 



Finally, the solution of (2.1 ) is given by the following table 



Ao < 



2o-2 



2^ 



Ao > 



2o-2 



a 



Ce""^^* 



+1 



a 









if 
if 



V'Ai.Aali) = ^ 
^Ai,A2(i) = ^ 
where a = \/\P' — 2X20'^ ja^ and C 

V'Ai.AalO 



> a. 



lAi -41 < a, 



^l^T" 



Ai-4y+a 



1-^(^1-^)* 






^Ai,A2(t) = ^ + /3tan (^/3^ + arctan (^^ 
where /3 = V2A2a^ - 67^2. 



2-b 
"2 — 



and t^ ^ is given explicitly in terms of Ai and A2 by 



10 





*Ai,A2 


1 ^^f Xj-b/cT^ + ^b^-2X2ayaA 
Vfe2-2A2CT2 \Xj-b/a'2-^b'2-2X2(72/a'2 J ' 


As < 1^ and Ai > ^ + ^/6^ 


- 2A2(t2/(J^ 


A2 < 1^ and Ai < ^ + V^ 


- 2A2a2/CT" 


t*\ \ = +00. 
Ai,A2 






^* _ 2 


"Ai,A2 o-2{Ai-^)' 


A2>|. 


*Ai,A2 ' 


2 r^ irctinr ^i'^"-^ ^^ 


V2A2<72-62 \^2 \yj2X2a^-b^ J J ' 



C Proof of Proposition |2.2 

/i!^(t) is defined as solution of 

t*x,M^l^\{t)) = t 

We distinguish 2 cases 

Case tl^,^,(M;) = ^,2,2'a2m;^^ ^"^ (xi^l^^T^^vSiST^) = ^"^ *^^' "^'"^ ^"^ ^^'''^' ^°' ° < 
e « 1, V'Ai(M;-e),A2K-.)(*) = ^-«£Szi-' ^here a = a(e) = ^^^ - 2A2(^; - 6)^2/^2 

and C = C{e) = ^ — r — ^ . We obtain, by writing tlie Taylor expansion of Ce'^" * 

with respect to e, that ip\Un* -A^\,jn* -e){t) can be written as 

b , .2-cie + 0{e^) 2a(0)/ci 
V'Ai(M;-.),A2K-e) W = ^ + «(0) cie + 0{e^) " ^ + '''^'^' 

where ei(e) is bounded. It follows that, for e small enough. 

On the other hand, we have 

</'AiK-e),A2K-e)(t) = j^* " "(e) (in (ce^* - e-^*) - In (C - 1)) (C.l) 

We obtain easily that V'Ai(/x* -e),A2(M* ~e)(^) = ^ In ^ + e2(e), where e2(e) is bounded 



11 



Case t 



Ai,A2 



V2A2a2-62 \ 2 



arctan 



^2A2<72-fe2 



. : In this case, for < e ^ 1, we have 



'/'Ai(/.;-e),A2{/.;-e)(0 = ^^ + ^ (ill COS ^(0, //;; - e) -In COS ff(t, //:;-€)) 



where 



g[t, /ij = 1 + arctan(- 



2 " ' V2^A2a2-62- 

We get by the same as before that for e small enough, 



where e(e) is bounded. 



e(e). 



(C.2) 



D Proof of Theorem [O 

D.l Upper Bound 

A simple application of Markov's inequality shows that, for every ^i > 0, 

F{Z > R)= P(e^^ > e^^) < e"'^^ E[e^^] = m e^+'^('^). 



In particular, for fi = fi* 



we have 



lnF{Z > R) < -in* - J^)R + A{fi* - J^) = -fj,*R + 2VojR + iyln\ - 

R y R \ uj 



It follows that, 



ln¥(Z> R) + fi*R-2VUR v 

'"^TJ ^ 2 



(D.l) 
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D.2 Lower Bound for (3.1) 



In order to derive a lower bound for the liminf for (3.1 ) we proceed as follows. Given any p > 
and any R > 0, we have 

= ^(j peP'-drlz<R]+E{eP^lz>R) 

< l + ^fy peP'-lo<r<z<RdA +E{ePHz>R) 

fR 

< 1 + / peP'^FiZ > r)dr + E (eP^lz>ij) 



(D.2) 



We choose p = ii* — y/uR ^, where /3 g] 2 , 1 [ to be given bellow. Denote by J := L e^'''¥{Zt > r)d. 



Using (D.l) we have 



R 






?l-/9|',/T_R'9-2\2 



Jo 

u+l^V^R^ ^-V^R^-^{l-R^-hf ^ ^gA(p)^z./2+lg-v^Ri-/5(l-H^-2)2 



Let's set 



2 \^/IjRlnR, 



For some positive 7]. We can then write R^ as 



rP = r\ e(/5-|)lnH ^ i^l (1 + (^ _ \ l^R) =r\{1+ h_^\ '- 

2 V vojR J 



It follows that (1 — R^ 2)2 is given as 



Thus 



We finally get, by taking ry = | + 4 



[l-R^'lf^l^ 



J < me^(P) R-2+^-'> 



E {eP^lz>R) > e^(P) (1 - R- 



(D.3) 



(D.4) 
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Applying Holder's inequality to E (e^ 12>_r), for some q > 1, we get 

E {eP^tz>R} < (IE e'lP^y^ (P(Z > i?))^"i . 
It follows that for any q > 1, we have 

(1 - -)lnP(Z > R)> InE (e*'^lz>i?) - -Hqp) 



Thus 



Q /^A/,„^ 1 A/„,„^ , 1.. /-, 7,-3^ 



We set 



lnP(Z>i?) > --^{A(p)--A(qp)+ln(l-R--^)] (D.5) 

g — 1 \ g 



j{R) = -^^, with 6g = -. (D.6) 



In particular, we have 



A{q{R)p) 



p ' '' R 

Hp) 



^, 



Thus 



'^(^) (A(,)-iAM) = ^#^fl-^)A(p) 



Hp)- 



q5q- 1 



5,{q - 1 
On the other hand, for R large enough, we have 



It follows that 

q5,-l ^ -^l^*{l - 5,) + (1 - 6l)V^R-^ - 6^{S - ly-jfi 
<?-! {1 - 6g)^R-l^ - 6{6 - lyj^ 

-fi* + (1 + Sg)^R-^ + 6^iy'j0- 
R-P + 5v^^ 
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Thus 



q6q- 1 



On the other hand, we have 



A{p) = -^* + (1 + 6g)Vi:jR-'' + 5'iy 



lni?^\ J^Rf^ + vlnRl^ 1 



R''^ J ^R-P + 5v^^^, 






1 ^/uJR'^ + u In Rf^ 

R2P 



Rf 



On the other hand, as -^ ~ i? ^ ^ we have -^ ln(l — R ^) = 0{R). We then get 



ln¥{Z > R)> {-fi* + {l + 5q)^R-^ + 6^v 



InRl^ 



) (^'"(^ - i?0 " ^^^^'^ 



Writing 



i?2/3 = ^ ^2,3-1 ^ ^ g(2/3-l)lniJ ^ ^ l + 2 



VuJR 



We conclude that for any a > |, we have 



l,mmf""''^^"'+"'^>0. 

R^ + OO i?" ~ 



(D.7) 



D.3 Lower Bound for (3.2) 



In order to derive a lower bound for the limsup for (|3.2|) we proceed as follows. Assume that 

u 3 



lnF(Z > R) + fi*R-2VujR 

limsup -— — = u" < 

R^oo HR) 2 4 



we will show that this will lead to a contradiction. 
We first write for p = n* — y^, 

e^(p) = E (eP^) = E f e^^l i "l + E f e^^l 



Z<R^ 



Z>R^ 



1 /-oo ^ 

pR^ I »^ / „-^/^z+2y/IJz+u'lnz 



< e^^' + m / ^ e"V ii 

'/?3 



dz 



Jr~^ 



where z/' e]i/*, | - 1[. 
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Lemma D.l. For any 7 G R, there exist c > so that 

poo 

/ ^ zT e-"^^^^''^^^ dz = cR-4 + 0{R-2) 
It follows that 



(Di 



which is impossible as z^' + | < | . 



Thus 



lnP(Z>i?) + /i*i?-2\VR z^ 3 
lim sup ; — --— > 7 . 

R^+oS HR) ~ 2 4 



(D.9) 



D.4 Proof of (|3^ 
Suppose 



hm ^ ^—4 = V* 

R^oo \n{R) 



we will show that v* = 
lead to a contradiction. 



i^ _ 3 

2 4- 



Assume z^* < | — |- For p = ji 



For this we will show that both {i^* < f - f } and {i^* > | - |} 
we have 



R' 



^^ gV^+f In(fi) ^ gA(p) < gpii. +M2 le' 



^z+2y^ujz+u' In 2^2 



jRr^ 



V' e-V^i^z-ifdz 



where u' €]i/*, | — |[. Using Lemma 



D.l 



this means that 



Mie^^i?^ < eP^'+cM2e^^i?"'+t~cM2e^^i?5i?-(5-|-') 



which is a contradiction. 

On the other hand, if we assume z^* > | — |, we will get 



Ml e^'^« i?f > e^^' + C2 M2 e^'^^ i?^'+i ~ C2 M2 e^'^^ fl? i?'^"-(f-t) 



withi/">G]zy*,|-|[n. 
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Proof of Lemma ID. II 

We note first that r°° 3 z"i e-v^^(^-i)'(iz = 2^3 z^"'^^ e'^^^^'-'^'^^dz. Denote by Ih) : 



'i?"3 



/J^S 



ji_3 ^7 e-v^^(v^-i)^^2 and 7(7) := J^ z^ e-'^^^^-^'>^dz. One can easily see tliat if 7 > 
/(H) > ^(7) > ^([7] + 1) and J([7]) < J(7) < Jib] + 1) 

and if 7 < 

/(H) < ^(7) < lib] + 1) and J([7]) > 7(7) > J([7] + 1) 

which means that for any 7 G M there exist ni, n2, ?^3, ?^4 G Z so that 



/(m) + J(n2) < /(7) + J(7) < /(ns) + J(n3) 



(D.IO) 



Now for any n G Z, we have 



/•I 

/(n)-I(n-l) = / ^ z^'-^z - 1) e-^^^'-'^^' dz 



R 



-5(^ „-V<:zr(/?-5-i)2 



/ujR 



VujR 
Using this equality we can easily check that 

/(n) =IiO) + OiR- 



n — 1 , 

+ ^=Iin-2) 
VZJR 



By the same way, we get 



where 



Jin) = JiO) + OiR- 






i-R- 



e-V^z^dz 



To calculate /(O) we consider the function loiR) := L e ^'^^^ d^;. This function is the 

unique solution of 



^o(^) 



[l + lRs)e 



^(l-i?"5)2yijR 



AR 



1 
AR 



loiR), /(1) = 0. 
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The solution of this equation is given as 



Io{R) = R-^ 



R-i 



1 /- {l + ly-l)e-('-y'^y'V^ 



Ay 4 



dy 



(l + iy-i)e-(i-2'"")'v^ 



4y4 



dy 



(l + iy-8)e 



-(l-y 5)2^/wj, 



. 3 

4y4 



where 



cq{I) R-'^ + 0{R~^) 



coin 



(l+iy-i)e-(i-?/"5)'v^ 



, 3 

4y4 



dy 



By the same way we get Jq{R) = ^Jo{R). Then 

JoiR) = co{J) R--^ 



where 



coiJ) = Jo(l) = / e 

'0 



We finally get for any 7 S M 
(co(/) + co(J))i?-3 + Oi{R-^2) < i(^) + j(7) < (co(/) + co(J))i?-3 + 02{R-"^: 



-dy 



(BAD 



(D.12) 



(D.13) 



E Proof of Proposition |4.1 



We have Fp{t) = Q"-Vp(t)+'"'4>p{t) ^ where ^p{t) = /^ iljp{s)ds and V'p is given by theorem 
Fp{t) = e."-fvit)+'"i'v(i)^ where ipp{t) = f^ipp{s)ds and V'p is given by 



2.1 



as 



2 ^ 2(72 



2 ~ 2cr2 



P^-P \ {b-papf 



> 



2^2 






V-plO 



b—pap 
„1 



O- IT — , II — V^\ ^ Ql 



where a = y^(6 — pap)'^ — p{p — l)a'^/a'^ and C 



Vp{^) - 2a^l+'-^Jf£t)- 



b— pap+acr^ 
b—pan—aa'^ 



Mi) = ^ + /3tan (/3^ + arctan (=^)) , 
where /3 = yjl^a^ - {h - papf /a"^ . 



where 7/;^ is defined for t G [0, t* [, where t* is given by 
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< p2 _ p < (J!_P^ and b - pap < 



p^ — p < - — '^2^' and b — pap > 



€ 



b— pap— W (b— pap)'^ —p{p—l)a'^ 
\/(b— po'p)^— p(p— 1)(t2 \b—pap+-\j(b—pcrp)'^—p{p—l)(T'^ 



In 



p^ — p = - — Sp^- and b — pap < 



2 ^ (b-po-p)^ 

P^ -p> 'I'Z ' 



t*p = +00. 



tX 



Ai,A2 b—pcjp' 



2(^1 



b—pap>0 



-arctan 



y p(p — l)cr-^ — (b — pqp)^ 
b — pap 



^p{p-l)u'^-{b-pap)^ 



Note that if p{p — 1) < 0, then t* = +oo. 

Consider the critical moment of e^. For t > 0, denote by 



^+(t) :=inf {p> : t* = t} and ^f (t) := - sup {p < : t* = t} 



(E.l) 



Denote by 
ti{p) := 



b — pap — yj{b — pap)"^ — p{p — l)a'^ 
^/{b — papY — p{p — l)o"^ V ^ ~ P'^P ~^ y^(6 — pap)"^ — p{p — l)a'^ 



:ln 



(E.2) 



and by 



t2{p) :-- 



2{.lb-p.p>o - arctan (^^^^^^Jf^^; 
Vp{p ~ 1)^^ — [b — papY 



i __ l-^±^(l-^)2+4(l-p^)^ ^^^^ ^^^^^^^_^ ^^ ^2 _ (6-p.p)^ 



Also denote by p^ := — - — 2(^"-n^^ ^ (^^^ solutions oi p^ — p 

have 



2(1-P^) 
.X 



//_ = - sup {p < Po : *2(p) = t} 



and 



/«^W 



inf {p > 1 : ti{p) = t} , If |p > 0, 6 - per < and ^r^— ^ < t\ . 



inf {p > Po^ : i2(p) = t\ , Otherwise. 



(E.3) 

). Then we 

(E.4) 

(E.5) 



Case /9 > 0, b — pa < and — 4—- < i : In this case, p^{t) is solution of ti(p^(t)) 



t. 



Which means that for e small enough and for any p G [/U+ (t) — e, /u:^'- (t)[, we have 



b - pap Ce"'^'* + 1 



a- 



a^ 



Ce^ 



j^t 



and 



W 



it) = ^-JP f - A (in (ce-^*/2 - e— ^*/2) - ln(C - 1 
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where 



2 



a = aip) = yih — papY — pip — 1)0"^ and C = Cip) = — ; ^ (E.6) 

—b + pap + au''' 



We can easily check that ij^x_^(t) and (p x_^{t) can be writen as 

i'.x.St) = ^ + 0(1) and (^^x_,(t) = ^ ln(i) + 0(1) 






where 

/ . 5 , — 1} r\i( 1 1 I 

General case : We have 

p^{t) := inf {p > : t2{p) = t} and n^it) := - sup{p < : t2{p) = t} 

where t2(j') := / 2 I In particular for any p G [P'j^{t) — €, p^{t)\, 

we have 

V'p(^) = n 1" /3 tan ( /3— — h arctan 



(T^ V ^ V per 

and 

h- pap 2 / arctan [-^^^ 



(Pp{t) = ^ — t + ^^ In , 

^ \/3^+ arctan (^ 



where 



/? = /3(p) = \l2^^a^ - (6 - pap)2/^' (Ei 



Denote by 

,(T^t , „ „x_ ( —h + pap 
~Ja' 

We can easily check that 



g{t,p) := /3^^ + arctan ( — ^-^ — ) (E.9) 



V^x_, = ^ + 0(1) and v.^x_, = ^ ln( J) + 0(1) 



where 
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By the same way, we show that for any p G] — /i_ , — ^_ + e], we have 

M^) = ~ r^ +/?(P) tan (g(t,p)) 



and 



It follows that 



where 



a^ 



b-pap 2 fcosg{0,p) 
fpit) = ^ — t + — In ' 



a^ 



0"^ \cos g{t,p) 



^-Mf- 



a; 



h 0(1) and w_„x 



(T 



^ln(Vo(l) 



W- := 






Note that uj± > 0, for any t > 0. 



(E.ll) 



F Proof of Theorem 14.2 



This theorem is very similar to Theorem 3.1 in |3] where having a sharp asymptotic formulas 
for the Call price as q-^i^+^2V^x'^, an asymptotic formulas for the implied volatility is obtained 
as in the theorem. More precisely, it is shown in [4] (cf Lemma 10.1 of ) that the solution a of 



-Aix+A2\/Xr^i^ 



X ~fc_).+c 



A^ 



-fc+|a?(fc)\ fc f-k-la^k) 



e^N 



Viat(k) ) \ yft(7t(k) 

utik) « h{A)^fk + /32(Ai)^2 + fi2{A^)iy + 



IS given as 

1 , log k 



2' Vk 



(F.l) 



where f3,{A) = ^ (^A+T-^) and (32{A) = ^^{^- ^). 



Vi 



Now suppose (3.3) holds for ibX(, then 

H > k) ~A;^+oo 



g-M+fc+2Vc^+x^^-4 and ¥{Xt < -k) ^k^+oo e 



-fi^\x\+2^/u)+\x\ 



a 3 



k^ 4 andEfe~^-e^l 



^fe— >+oo 



In particular, we have E (e — e ) ~a,-->+oo e 

e~^ '''^- ' ^ A;o^~4 . It follows that at{k) and at{—k), for /c large enough, are given as solution 



to e-(^+-^)'=+V<'^fc^-i ;^ AT / --S-. 



,-feAr ( -k+li<^ti-k) 



e-^N 



VitJt{-k) 



Vt(Tt(k) 

N ( -fc-l^t(-fe) ^ , Thus 



Vto-t(fc) 



, a 1 , log k 



at{k) ^ Piifi^ - l)Vk + 2/32(/x:^ - 1)V^ + P2{f^t - l)(o - t) 



2 i' Vk 
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and 



o-i 



(-fc) « /3i(A.^)V^ + 2/32(/i^)V^ + /32(Ai^)(" - -) 



^2 A' ^ 

Now as the Call price is an increasing function of the implied volatility, an upper bound of 
the distribution given as g'^'^+^v'^^x''* will give an upper bound of the implied volatility. So we 



can easily prove (4.4) this way. 
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